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ABSTRACT 


Utilizing four-quadrant experimental data determined by Minio-~ 
Vich converted to the modified advance and angle-of-advance 
coefficient form, complete Karman-Knapp Circle Diagrams are 
plotted for a number of 3 bladed propellers. Reducing’ the 
family of curves to a rational basis, algebraic equations are 
derived to model the performance of a three bliaded propeller 
through ali four quadrants using Eulerian turbomachine theory. 
The result is a set of equations which can be used to generate 
complete propeller thrust and torque characteristics given only 
incomplete performance data suitable for use in steady state 
design or dynamic simulation. 


Thesis Supervisor: Dr. Henry M. Paynter 


Titles: Professor of Mechanical Engineering 
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CHAPTER 1 INTRODUCTION 


1.1 Primary motivation for this Kind of modeling method 


For a naval architect or anyone interested in the modeling of 
Propellers as a part of the design procedure leading to an 
integrated ship design, the problem long persisted where one 
1s required to estimate the performance of a propeller’ for 
ship modeling purposes before the design has progressed to the 
point where one 1s reasonably confident as to what the final 
desired propeller should look like. Even more frustrating is 
the problem of having somehow determined the geometry of the 
propeller but not having sufficient experimental data to allow 
a comprehensive computer simulation. 

Even where it has existed, considerable effort has been expen-— 
ded to convert the data into useable form for computer simu- 
lations since the practice of tabulating the results as func- 
tions of an advance coefficient has difficulty in the vicinity 
of zero ship or shaft speed. 

The problems associated with modeling a machine such as a 
propeller are formidable since a complete model must make 
allowances for complicated fluid motions which are only just 
beginning to be understood. In particular, it is still quite 
difficult to pass completely from propeller geometry to the 
detailed lift and drag forces required for specific propeller 
characteristics. 

However, there now exists a considerable amount of documen- 


tation for modeling machines very much like a propeller. For 





example, industrial axial turbines and pumps have had _ the 
attention of modelers and those who perform simulations for 
decades. TaKen together, propellers, turbines and pumps are in 
fact, members of a general class of machines Known as turbo- 
machines. Models for many of the Jarge industrial machines 
have been quite successful in predicting machine performance 
characteristics and are widely used in control system = synth- 
esis and design. A general model for propellers which could 
incorporate that data base would have a firmer basis in exper-— 


ience and first principles. 


i.2 Some current propeller modeling methods 


At present there are two common methods used by designers’ to 
obtain propeller data for computer simulations. The first 
involves experimental determination of the characteristics of 
a physical model of the full size propeller followed by con- 
version of the data into the familar thrust and torque coeffi- 
cient form. Such coefficients are then used in a table look-up 
scheme within the simulation program to obtain the necessary 
thrust and torque values. There are advantages to this Kind of 
method. First, the determination of the propeller characteris-— 
tics through experimentation (scaling considerations aside) 
ensure that the data used in the simulation will be as accu- 
rate as possible. Secondly, the use of a model of the intended 
propeller acts asa test of the Be onneey and reveals at an 
early date any serious conceptual problems. 


The disadvantages of such a method are fairly obvious. To test 





a model of a propeller at an early stage in a design, you must 
have a well defined idea as to the range of operation of the 
propeller and its desired geometry. To allow for a comprehen- 
sive simulation, you must test and record the model] data in 
all four quadrants of operation. Since cost is always a signi- 
ficant consideration, the machining of a propeller model] and 
the set-up and operational costs of the tests represent a 
sizeable investment in the design. As alluded to earlier, raw 
data is not suitable for computer simulations due to the 
problem of ‘singularities’ of the advance coefficient at zero 
ship speed and zero shaft rpm. The experimentally determined 
data must be converted to a form that avoids such difficul- 
ties. 

The second method of modeling propeller characteristics invol- 
ves complex and relatively recent work in propeller hydrody- 
mamics. The computer models of Professors Kerwin [1] and 
Cummings [2] are quite thorough and provide some revealing 
results. As a ’complete’ model, these computer programs are 
undoubtedly at the vangard in propeller theory. With a firm 
Basis in hydrodynamic theory, the continual evolution of these 
programs will ultimately achieve very accurate results. Unfor- 
tunately, these programs are immensely long and require consi- 
derable amounts of CPU time. Because of the attendant expense, 
they are not yet practical for incorporation directly into a 


ship simulation program. 





\ 


1.3 The Objective of the Modeling Equations 


The objective of this thesis will be to derive a set of rela- 
tively simple algebraic equations that will serve as a_ tool 
for naval architects in the ship design process. By being able 
to predict accurately propeller performance in all four quad- 
rants of operation based on Knowledge obtained from !imited 
experiments and iin a form suitable for computer simulation, 
the naval architect can save considerable time and expense in 
the simulation portion of the design. Experimental determina- 
tion of the points of minimum dissipation would be sufficient 
to allow for a comprehensive simulation for forward or crash- 
back maneuvers. Such a set of equations may not eliminate the 
need for actual experiment in the final stages of the design, 
however they would greatly simplify the initial requirements. 
Additionally, where experimental data exists for particular 
propellers in only the first quadrant, the equations could be 
used to reconstruct the remaining three quadrants with good 


accuracy. 





CHAPTER 2 DERIVATION of the ALGEBRAIC MODEL 


2.1 Fundamentals of Eulerian turbomachines 
The practice of presenting open water propeller performance 
data as nondimensional thrust (Kt) and torque ¢Kq #£coeffi- 
cients versus an absolute advance coefficient (lambda) = or 
simply advance coefficient ¢€J> (31,f4] incorporates assump- 
tions which inherently restrict the machines of interest to 
Eulerian turbomachines. That 1s, by considering only high 
Reynolds numbers the performance of the machine is governed by 
inertia and pressure forces. Dimensional analysis on the ex- 
pression: 

(Thrust) F=#¢(Vp,D,N, P> 
yields 

Kt= F/PDINZ=F¢(Vp/ND) =F¢ UN) (2.1) 


Similarly for torque, 


Kqg= T/PDENZ2=Qg(Vp/NDD =Q6 Ud (2.2) 


where Vp is the speed of advance of the propeller through the 
water and D any characteristic propeller dimension Cusually 
diameter). This single parameter representation of propeller 
performance yields the powerful result that the performance 
data of geometrically similar propellers will, when plotted 
together, generate to good approximation a single curve. 

The same assumptions would , from a fluid mechanics stand- 


point, give Euler’s equation of inviscid motion. That is, 


PYEGRADCY) = —-GRADCP) 2.3) 


10 





Or, in a more general form to incorporate unsteady flows, 


PDOCVISDt = -GRADCP) (2.49) 


where D/Dt is the substantive derivative. 
In conjunction with equation (2.94), we express mass conserva- 


tion by: 


DCP) /Dt + PGRADKY = 0 (2.5) 


where P is independent of position andv is the vector fluid 
velocity. Reference [5] chapter 3 and Reference [6] are useful] 
for amore detailed analysis. 

A machine satisfying equations (2.94) and ¢€2.5) will, as veri- 
fied experimentally, demonstrate what is referred to as 
"Eulerian similitude". That is, geometrically similar machines 
will produce similar performance characteristics when present- 
ed in a properly chosen nondimensional form. Casanady [7] 
refers to this as a “performance law’. 

Based on the discussion above, the determination of the per- 
formance characteristics of a large turbomachine is usually 
done with small scale models. The full scale results are 
derived from the model results by recognizing that at corres- 
Ponding points on a characteristic curve, all nondimensional 
Parameter combinations are also equal. 


However, there are problems with scale effects not evident in 


ed 





such a single parameter representation. In addition to geo- 
metric dissimilarities, such as roughness and clearances, hy- 
drodynamic differences are also present. These include Rey- 
nold’s and Froude number variations as wel] as cavitation and 
compressibility number differences. Attempting to incorporate 
some or all of these variations into a performance law would 
require that the law be multidimensional. As an example, one 
could attempt to formulate a nondimensional representation for 


thrust and torque characteristics by: 
Kt=#¢J,Red Kq=f¢J,Re> . 


Here, the equivalence of the advance coefficient in conjun- 
ction with the Reynold’s number would lead to similitude. 
Similarly, one could also form yet another expression for 
thrust and torque 
Kt=f¢J,P) Kq=f¢J,P) 

where P represents a nondimensional] ratio for propeller pitch. 
Such a formulation could be applicable to Kaplan turbines’ or 
controllable pitch propellers where the blade pitch is varied 
mechanically. 

Multi-dimensional performance laws unfortunately lead to fami- 
lies of SUIpECe. a problem which single parameter laws avoid. 
AS an example, Bovet [8] presents efficiency versus output at 
a constant head for a Kaplan turbine for various blade set- 
tings. An analytical treatment for a multi-dimensional per- 
formance Jaw in terms of generating numerical data would 


entail a matrix of possible states with dimension equal to the 


He 





number of nondimensional characteristic parameters. Obviously, 
dimensional analysis techniques used to peneate characteris-— 
tic nondimensional performance parameters are limited in their 
applicability in performance simulation. To allow incorpora- 
tion of fluid or geometrical affects into the analysis would 
require another approach. That approach should produce a for- 
mulation which retains the advantages of similitude yet allows 


for some fluid or geometrical effects. 


\ 


2.2 Ideal, canonical and perfect ship propeller equations 


Figure 2.1 shows a propeller pics cross section at some 
arbitrary radial distance from the propeller hub under’ the 
condition of zero nominal speed of advance. Figure 2.1la in- 
dicates the tangential velocity of the blade element (Vr) and 
fluid velocity (Vf). Figure 2.i1b shows the lift ¢(L2 and drag 
(D) forces with respect to an effective angle of attack a The 
magnitude of the effective angle of attack is a function of 
the circulation which causes an induced velocity [3]. 


If we define a torque exerted by the element as Fr then, 
Fr= LSint@m + DCost(m@ (2.6) 


Further, defining Ff as a force perpendicular to Fr in the 


direction of fluid motion (fluid pressure term, 


Ff= LCos(m - DSin(@ (2.7) 


13 








Ca? 


CB) 


Figuge 2. 1 
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Since L and D are expressed as: 


L= C1C%PAVEY) 3 D= Cd ¥PAvxy) —, 


(2.6) and (2.7) can be rewritten: 


Fr= 4PAVCCIEVE + Cdk¥r) (2.8) 


Ff= 4PAVCCdkVr — CIRVF) (2.97) 


In per unit (see Appendix AD, (2.8) and (2.9) become: 


m= V(C1¥q + Cd¥n) (2.10) 


p= vCClkn - CdkqD (2.11) 


Or using the canonical notation, 


m= RgxXq + Rdk&n Cet?) 
p= Roxkn —- Rdk&q C2. 13) 
where 
Rg=r qkv rq=rq¢ F) V=SQRT (nkn+qxkq> 
Rd=r dkv rd=rd¢F) P=ARCTAN( q/n) 


The correspondence with the bond graph representation shown in 
Appendix B is now readily apparent. 
The mechanical power can be expressed as Frx¥¥r while the fluid 


power can be expressed as FFXVF. Thus, dissipation losses can 


1S 





be found to be 


Pd= Frx¥Vr — FxVF €2.149) 


In per unit, ¢€2.14) becomes 


Pd= mkn - pxkq (2.12) 


Substituting ¢€2.12) and (2.13) yields 


Pd= Rdiviv €2.16) 


Hence, as would be expected, an ideal ship propeller would 


have no power lost in dissipation if it could generate lift 


without drag. The equations reduce to the rather trivial 


expressions 


J 
il 


Rgxq €2.1/7a) 


RoXn (2.17b) 


vv 
bi 


with Rq representing an energy transformation modulus. Rd can 
be considered as a general loss term that includes a combin- 
ation of a number of loss mechanisms such as form drag and 
separation losses. Writing Rd in a more general form typical 
of such losses such that at the point of minimum dissipation 
loss, Rd is identically zero [9], we can define a PERFECT 


propeller. 
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Rd = WOR) ¥CCn-gmk(n-gm/yv) - (2.18) 


Using the Kutta condition to formulate a general expression 


for the modulus Rg would suggest 


Rg = rq + XC PH) ¥CCn¥nd -Cqkq)/v C2017) 


differing from the form proposed by Paynter By the term rg 
which could represent some linear function inn. Substituting 
(2.18) and (2.19) into (2.12) and €2.13) and recognizing the 
arbitrariness in WCF) and XCF), we can arrive at the following 


expressions: 


3 
ll 


rQgoxq + K1i¢Ch) xv¢n-q) €2.20) 


Pp = rgkn + K20F) ¥vCn—-q CZn 212 


where Ki and K2 now incorporate the individual character of a 
given propeller. The equations above for a PERFECT propeller 
(zero minimum dissipation loss) can be augmented to represent 
a PRACTICAL propeller by including an additional term to 


reflect the actual non-zero minimum dissipation loss. 


m= rgkiq + KiCPHkvtn-q + RiXn €2.22) 


rgkin + K20F)kVvCn-q - R2Xqg €2.23) 


U 
" 


Ri and R2 would be functions of the advance coefficient J for 
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a given propeller, but in a simpler model could be reduced to 
constants. As will be seen later, the values of RI and R2z have 
a direct effect on the shape of the 4-quadrant curves of a 
Given machine and in particular, the slope of the zero torque 
and thrust radials. K1¢® and K2¢F) for an idealized propeller 
would be identical. The degree of agreement between Ki and K2 
based on experimental data would be a measure of the validity 


of the model for a given rq. 


2.3 Determination of the alqebraic modeling equations 


Using equations (2.22) and €2.23> as a model for determining 
propeller characteristics requires a determination of the K¢CP) 
term as well as R1, R2 and rg based on experimental data. The 
data produced by Miniovich [10] and converted to a modified 
advance coefficient form [11] was used in an abbreviated 
methodology in order to formulate a general model that could 
be expanded upon in amore detailed analysis. See also Refer- 
ence £12]. Such an analysis would include an investigation of 
those factors which collectively yield the individual charac- 
ter of a given machine yet have defineable specific influen- 
ces, the details of which are beyond the scope of this thesis. 
In Reference [C11], experimental results for thrust and torque 
are tabulated for eighteen three bladed propellers with vary- 
ing expanded area ratios and pitch-to-diameter ratios. The 
propeller designations are those used by Baker and Patterson. 
Table 2.1 is provided as a quick reference for EAR and P/D 


values for each propeller. 





Figures 83 and B4 in Appendix B are plots of raw data for two 
series of propellers. To be useful, these curves must be 
rationalized to a consistent basis. To do so, the points of 
minimum dissipation were found for each of the propellers (See 
Appendix A). Those points were then used to scale figures B3 
and B4 such that the point of minimum dissipation in the 
region of normal pump operation (quadrant I> was fixed at 
€1,1> and in the region of normal turbine operation (quadrant 
III) at C-1,-1) in the qvs n plane. The mapping of the 
data to those points of minimum dissipation is done over’ the 
entire qvs n plane recognizing the “mirror image’ pattern of 
the curves. That is, the point of minimum dissipation in the 
first quadrant for positive torque and thrust is reflected 
into the third quadrant on the negative torque and thrust 
curves. Hence the scaling is reflected through the origin 
rather than simply across the zero torque and thrust radials. 

Figures 2.2 and 2.3 are plots of the thrust and torque curves 
respectively for three propellers selected for their varying 
expanded area and pitch-to-diameter ratios and are representa-— 
tive of the form of the 4-quadrant curves for all of the 
propellers. With the data plotted in this form, the similarity 
of the torque and thrust curves versus the advance angle is 
striking. In addition, the general form of the curves is quite 
consistent in the regularity of the specific perturbations 


which influence all of the propellers in the same manner. 


i 








Figure 2.2 
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Prop No. 1 2 3 a 2 6 
EAR GS Geo "G.5 0.5 6.5 0.5 


P/D 0.6 0.8 1.0 1.2 1.4 1.6 
Prop No. 7 8 ? 10 ii 12 
EAR 0.8 0.8 0.8 0.8 0.8 0.8 


P/D 0.6 0.8 1.0 1.2 1.4 1.6 


Prop No. 13 14 i5 14 17 18 
EAR cote. ee. 1 


P/D 0.6 0.8 1.0 1.2 1.4 1.64 


Table 2.1 


Before continuing, I must point out that there 1s an 
approximation which was used which requires an explanation. As 
Figure B2 in Appendix B demonstrates, the lines of zero 
thrust Chead) in the first and second quadrants are obviously 
not 180 degrees apart. The difference between the angle 
measured clockwise from the first quadrant zero thrust radial 
and the second minus 180 degrees I call the sKew angle. The 
skew angle for high specific speed axial flow pumps is 
generally quite small. See reference [13]. The skew angle 
computed for the propellers used in this analysis were all 
less than 3 degrees and were generally less than 1 degree. The 
amount of skew angle is a function of the magnitude of losses 
of the machine. Since the sKew angle is quite small, the 
curves were approximated to a PERFECT propeller by dividing 


the skew angle between the two radials found from the minimum 
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dissipation calculations. 


then 
the first 
apart. The 


lated minimum dissipation angles were obtained with a 
dence of plus or minus 2 degrees due to the required extrapol- 


ations in the corresponding portion of the experimental 


Such an 


modeling 


calculated minimum dissipation angles 
angles in 


along with the calculated and approximated magnitudes of 


thrust 


third quadrant asymptotes were 


equations 


considerably. 


and torque radials at those 


values are shown in parenthesis. 


Prop. No. 


13 


18 


Min. 


Dissipation 


angles 


35.2 


39.5 
-142.5 


44.5 
-136.0 
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-134.5 
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44.3) 
135.8) 


33.0) 
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45.46) 
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recalculated using these adjusted angles. 


approximation also simplifies the derivation of 


Table 


and 


both quadrants for six of the eighteen 


2.2 


the 


The scaling values for q and n were 


As a 


180 


approximation was felt justified since the calcu- 


lists 


angles. The approximate 
Radial Length 
thrust torque 
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7.-498(7.07) 20.99¢€18.76) 
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Using equations (2.20) and (2.21) for the PERFECT propeller 


and setting rq equal ton, we can plot the KCF) functions 
required to fit the experimental data by setting m and p to a 
constant ¢1) and solving for Ki and K2 at each data point. 
Since Ki and K2 must be positive at all times, we account for 
the sign change inm and p in the different quadrants as the 
calculations proceed. By plotting Ki and -K2 as functions of 
the normalized advance angle, we can gain a quick qualitative 
feel for the similarity of Ki and K2. Figures 2.94, 2.5, and 
2.6 are plots of Ki and -K2 for propellers %?, i135 and 18 
respectively. 

As is readily apparent, even with a simplified function for roa 
(nd, Ki and K2 are very similar. The behavior of the KCF) 
functions in response to changes in pitch-to-diameter and 
expanded area ratios is also apparent and are manifest in the 
amplitude variations and relative smoothness of the plots. 

On the basis of that similarity, we form a new pair of KCF) 
functions by the straightforward transformation: 


Kal F (Ki + K2)/2 


KB¢ PD 


CKi - K2)/2 

such that Ki equals Kat+Kb and K2 equals Ka-Kb. Hence, Ka 
represents an average of the two functions while Kb is a 
measure of their difference. Figure 2.7 is a plot of the Ka 
and Kb functions for propeller number 15. 

A different approach may be to account for the general ampli- 
tude differences in Ki and K2 on eitherside of the zero thrust 


and torque radials and themselves. Four gain parameters could 
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be obtained for each propeller such that when Ka and Kb are 
plotted, the maximum amplitude of KBCF) will Be considerably 
reduced. This was done for propeller number 15 by gain scaling 
the average amplitudes as outlined above. Figure 2.8 shows the 
result for Ka and Kb. KbC®) now approaches what may be con- 
sidered as individual propeller nuance characteristics. 

Since Kat is essentially doubly periodic over the range of F 
we form yet another function that can be used to approximate 
Ki and K2 valid over the entire range of F . By overlaying the 
plots of Kat for -135 ¢< F< 45 and 45 < F< 225 and 
again averaging we obtain an approximate K(®) function applic- 
able to both the thrust and torque curves and valid for the 
two ranges of F stated above. 

To retain all of the information available such that we are 
able to reconstruct the original Ki and K2 functions as 
necessary, we must define four additional functions. Those 
functions are formed by first separating Ka and Kb on the 


intervals of -135 < F< 45 degrees and 45 ¢< F< 225 degrees. 


Labled: 
Kal = Ka -135 < F< 45 degrees 
Ka2 = Ka 45 < F< 225 degrees 
Kbil = Kb -135 < F< 45 degrees 
Kb2 = Kb 95 < F< 225 degrees. 


We then define the four functions that represent the average 


and differences of these functions on the two intervals of BP. 
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K = (Kal+Ka2)/2 
1 
K = (Kal-Ka2)/2 
Zz 
K = (Kb1it+tKb2/2 
3 
K = (Kb1i-Kb2)/2 
4 
Algebraic manipulation will show that: 
Ki = K +K +K +K -135 < F< 45 deg 
123 4 
Ki = K -K +K -K 45 < F< 225 deg 
123 4 
K2 = K +K -K -K -135 < F< 45 deg 
123 4 
K2 = K -K -K +K 45 < B < 225 deg. 
123 4 


Our approximate KC) function is now found from: 

KC ® = (CK +#K 0/2 

i 3 

valid on the intervals of -135<FK45 and 45<RK225 cegrees. 
If the two functions K1 and K2 were the same and the identical 
underlying function was responsible for the thrust and torque 
characteristics on each side of the 4-quadrant asymptote, then 
K , K , and K would be zero. In summary, KC® is an average of 
ce a K2 Nese is then symmetrized for the + and - torque and 
thrust curves. Figure 2.9 shows the result for propeller 
number Y on the interval -PI < F< PI. 
To determine some of the underlying fundamental characteris- 
tics of the KC function we can formulate an expansion of 
KCF). Two of the methods which are available involve the Walsh 
functions and the Fourier series. With the Walsh functions, 
square waves are superimposed and in the limit can be used to 
approximate simple sinusoids. They could be particularly use- 


ful here since the KC® function does have very sharp minimums 


at the $5 and -135 degree angles. The Fourier series is some- 
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what more familar and was used here. As a dual to the Walsh 
functions, in the limit, the Fouries series can be used to 
approximate square waves. 

By obtaining the Fourier series expansion for KC) we have 
finally, a potential model for this family of propellers. 
Table 2.3 lists the first five Fourier coefficients for 
propellers 9, 15 and 18. The Fourier expansion is taken over 
the interval between the minimum dissipation angles. At those 
angles the KC® function approaches zero and provides a 
convenient point to begin the expansion. The expansion i185 
Qiven by equation (2.24 where c 1s defined by equation 
(2.25) . 

By substituting the K¢® expression into the basic modeling 
equations €2.20) and (2.21), we can compare the experimentally 
determined thrust and torque curves to the model by plotting 
the equations over the range of -1 ¢<¢ n ¢ +1 . #=+%IIn doing so, we 
are assuming that the differences in Ki and K2 are sufficient- 
ly small to permit their omission. Justification of this 
assumption lies in the accuracy of the curves generated by the 
model equations when compared to actual experimental data. 
Figures 2.10 and 2.11 show the results for propeller number 
is. 


KC F) aQ + SUMC atk) Cost(Kkx) + BCKOSinCKx) )D €2.24) 


SQRTC aka + bb ) €2.25) 
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Propeller Number 


9 1S 18 
a 5 Cc a b Cc a b Cc 
k=0 ».20 0 3.20 6.27 0 6.27 10.31 0 10.31 
kK=1 o0O7 2.26 - 64 64 Say 87 ae £95 re 


Pa2 —-2.02 -1.06 2.28 —-2.68 -1.15 2.92 -3.27 -1.2? 3.92 


K=3 -1.25 2/7 1.28 -~.07 -~-19 -60 —-1.25 12 1.26 
K=4 -~.9? -.87 1.06 /3 -.223 24 -.78 01 7s 
Table 2.3 


To demonstrate the power of the model the curve resulting from 
utilizing only the zeroth term of the Fourier series expansion 
is plotted along with the curve plotted using the first three 
terms. Doing so reveals that the difficulty in fitting the 
proper K function lies essentially in the dissipation regions 
of the 4-quadrant curves. In a sense, this is expected since 
the KC™) function represents a modulating dissipation term in 
the basic modeling equations. What is interesting is that the 
plots indicate that a common dissipation mechanism affects 
both the torque and thrust characteristics and are lumped 
together in the K¢C® term. 

Additionally, it is obvious that quite accurate results for 
normal ahead propeller operation and for turbine operation can 
be achieved with K¢F) assumed to be a constant. As more terms 
of the expansion are added, the fit can be improved in the 
dissipation regions. For a loss of power simulation for exam- 


Ple, the most accurate results for thrust and torque values 
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would require a large number of terms. The limitation on the 
number of terms that can properly be added-to the expansion 
lies in the residue of the Ki and K2 functions. Recalling that 
KC) was formed from an averaging process involving Ki and K2, 
the functions which represent their differences were formed to 
allow retrieval of the original Ki and K2 data. To the extent 
that these differences degrade the KCP) = approximation, their 
impact should properly be less than the last term added to the 
expansion. AQain, the Walsh functions may reduce the number of 
terms required particularly in the normal pump and turbine 
areas of the 4-quadrant diagram. The model equations can also 
be used to construct the familar Kt and Kq versus J curves. Kt 
is given by equation (2.26) and Kq by equation (2.27). J is 
simply q/n where R= nkn+qkq and THETA=Arctan J. 


Z 


Kt 1./RXcos (THETA) (2.26) 


2 


Kq 1./Rxcos (THETA) <2uct, 


Here, the R andn refer to the appropriate values for the 
thrust and torque curves. 

By selecting ann and solving the modeling equations (2.20) 
and ¢€2.21) for q with a given KCP), Kt and Kq can be quickly 
determined for any value of J. Figure 2.12 is a plot of the Kt 
and Kq curves for propeller number 18 using the zeroth term of 
the KCF) expansion and the first three terms for J > zero. 


Note that the scale factors used in constructing the rationa- 
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lized plots of thrust and torque must be considered in calcu- 
lating Kt and Kq so that their ranges of values are those 
normally encountered in the literature. The scale factors used 
for propeller number 18 in constructing Firgure 2.12 can be 


deduced from Table 2.2. 
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CHAPTER 3 INTRODUCTION of ADDITIONAL TERMS and OTHER 
DETAILS in PROPELLER CHARACTERISTICS 


3.1 Pitch as a primary determinant 


As will be discussed in Chapter 4, the variation in pitch 
expressed as a pitch-to-diameter ratio (P/D) has a significant 
effect on the performance characteristics of a given propel- 
ler. From ae designer’s standpoint, Saunders [3] addresses 
pitch-to-diameter selection primarily on the basis of obtain- 
ing maximum propeller efficiency, particularly when utilizing 
series chart data in the initial stages of the design. From 
the perspective of the modeling equations however, it 1S more 
convenient to see the relationship between the relative pitch 
to diameter ratio and specific speed. References [14] and [15] 
present descriptions on the use of the specific speed para- 
meter. 

Specific speed variations have a significant effect on the 4- 
quadrant diagram not easily seen in the propellers presented 
here. Since propellers are in general very high specific speed 
machines, the impact of variations of that parameter is some- 
what obscure. The principle effect of specific speed variat-— 
ions on the 4-quadrant diagram was only touched on in the 
plots of the thrust and torque characteristics discussed ear- 
lier. It 1S seen primarily in the location of the zero thrust 
and torque radials of Table 2.2. As the relative specific 
speed increases, the radials tend to rotate counterclockwise 


while the skew angle addressed in Chapter 2 tends to decrease. 
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How variations in pitch affect the specific characteristics of 
a propeller is best seen in an analysis of the performance of 
controllable, reversible pitch propellers (CRP). Reference 
[16] is quite thorough in examining quantitatively the impact 
of pitch on the efficiency and thrust and torque coefficients 
of a CRP propeller. It is seen in the modeling equations as a 
dominant factor in the determination of the amplitude of the 
KC) function. An increase in the pitch-to-diameter ratio has 
two noticeable effects. First, it tends to increase the rela- 
tive amplitudes of the various harmonic components in K(CPF) 
which results in a somewhat wider fluctuation of the 4- 
Quadrant curve in the principle dissipation regions. Secondly, 
it has a direct affect on the zeroth term of the Fourier 
expansion which represents the average value of the K¢® fun- 
ction. The net effect is a shrinking of the 4-quadrant curves 
Ccorresponding to larger torque and thrust coefficients). 
Obviously, the “dc” value of the expansion is dominant since 
an examination of Figures 2.2 and 2.3 would indicate that the 
harmonic content of K¢(F) in the dissipation regions, specifi- 
cally in quadrant IV is apparently proportionally equal for 


all three propellers. 


3.2 Number of blades and blade area 


The effect of variations in numbers of blades was not examined 
in detail. However, as pointed out in Chapter 4, the formula- 


tion of the modeling equations was quite general and utilizes 
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what 1s recognized as the actuator disk theory to obtain their 
initial relationships. A comparison of 4-quadrant diagrams for 
three bladed and seven bladed propellers has not indicated any 
potential problems. What is really the primary determinant is 
the amount of blade area expressed here as expanded area 
ratio. Yet Figures B4¢a), B4¢b), 2.2 and 2.3 would seem even 
to discount any major influence of EAR on torque and thrust 


characteristics in general. 


3.3 Effects of cavitation | 


Cavitation 1s probably the most difficult aspect of propeller 
operation with which to deal. Rightly so, a lot of time and 
energy is being expended to describe the mechanics of cavita- 
tion and its influence on propeller characteristics. Referen- 
ces (4], [16], £17], €18), and £19] explore qualitatively and 
quantitatively the role of cavitation in propeller perform- 
ance. Reference [16] relates the speed of advance of a propel- 
ler to the cavitation number. By defining the cavitation 


number (sigma) with equation (3.1), 


sigqma=2qgChthA-hV) “VakVa CS. 12 


where g is the acceleration due to gravity, h the depth of the 
propeller axis beneath the water surface, hA the atmospheric 
head pressure, hV the vapor pressure of water and Va_ the 


Propeller speed of advance, it is shown that Cusing typical 
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values) sioma is related to Va by the relationship of equation 


(3.2). 


signma= 1070/VakVa (3.2) 


with Va in Knots. How the cavitation number affects Kt and Kq 
is seen in Reference [4], Chapter 16, Figure 108 as a fall off 
of Kt and efficiency as the advance ratio decreases. It does 
appear that the effect of changes in sigma for a higher J is 
smaller, also seen in Reference [16]. 

Yet nowhere is the role of cavitation included directly in the 
modeling equations. Since the intention of the modeling equa- 
tions 1s to provide a reasonable method of determining propel- 
ler characteristics at an early stage of a ship design say for 
the purposes of simulation, perhaps an extension to include 
cavitation affects is unwarranted. However, we can anticipate 
where such affects may be present. 

Equations (2.22) and (2.23) can be modified to indicate gen- 


eral relationship to cavitation number shown below. 


m=rg(sigmadikgq + K1ICF,siqmad kvk(n-q + R1i¢€sigqma) kn (3.3) 


p=rg(sigmad)kn + K2¢0F,sigma) ¥vk¥(n-q) - R2¢sigqma)kq (3.4) 


In some cases, the presence of cavitation is handled by an 
effective reduction of the propeller blade area. This could be 
incorporated into the K(f =functions. However, the greater 


dissipation would alter the 4-quadrant zero thrust radials to 
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some extent, hence, the Ri and R2 terms are included in this 
case. Notice that R1 is also a function of sigma Since the 
losses will be reflected from the hydraulic to mechanial side. 
Although the potential is there for inclusion of cavitation 
affects, such a step may be beyond the scope of the assumed 
model equations. Again, the intent of the equations is to 
provide a relatively simple means to simulate propeller per- 
formance for first-cut design purposes. Any extension beyond 


that goal will increase the complexity of the equations rapid- 


ly. 


3.4 Hull/propeller interactions 


Two of the familar parameters which reflect the influence of 
hull/propeller interactions are the thrust deduction factor 
Ct) and the wake fraction (wd). The wake fraction appears as a 
consequence of the disparity in ship forward velocity and the 
velocity as seen by the propeller. The wake fraction is de- 


fined as: 


w= (V-Va) /V 


where V is the ship velocity and Va the velocity that the 
propeller sees. 

The thrust deduction factor results from the requirement of 
the propeller to produce enough thrust to propel the ship 


through a resistance of Rt while also reducing the high pres- 
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sure area over the ship’s stern. Defined as: 


t = CT-RtO/T 


where T is the net required thrust to propel the ship at a 
given speed. Values for both factors are usually obtained from 
model tests. Reference [16] would indicate that t and w vary 
little over a wide range of ship speeds. In terms of the 
modeling equations, the approach was one of considering the 
propeller in an open water condition and hence, corrections 
for t and w were not included or necessary. Owing to the 
actuator disk approach essentially taken here, allowances for 
complicated inflow fluid geometries due to disturbances ahead 
of the propeller have been ignored. 

Reference [16] also addresses the pitch dependency of the 
thrust deduction factor and gives general ranges of t for P/D 
ratios of ti to -1i. References [4] and [20] address the incor- 
poration of the thrust deduction and wake fractions in more 


detail. 


3.5 Other considerations 


AS formulated, the modeling equatins simply do not address 
many of the fine details of propeller design. Most of these 
however, including blade rake, skew-back, the shaping of blade 
edges, hub fairing and blade profiles are in large measure 


details which go beyond the requirements for the initial 
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stages of propeller design. Once the ’macroscopic’ character- 
isics of a generic propeller type is nate naieec the impact of 
cavitation and specific blade loading requirements can be 
examined but only in the detail allowed by actual model tests 
or by a complex computer program specializing in propeller 
hydrodynamics. 

A tacit assumption which has been incorporated from the begin- 
ning is that the nondimensional, per unit quantity gq is ob- 
tainable from the dimensional values of propeller speed of 
advance Va and in turn, the flow rate through the propeller. 
There 1s, however, a term called the axial inflow factor a [4] 
which 1s included in the momentum theory for propeller action 
not directly accounted for in the calculations. It has been 
assumed to be a quantifiable constant. 

The efficiency of a given propeller and its relationship to 
the behavior of the KCF) function has not been explored. The 
calculation of the efficiency of a given machine is straight-—- 
forward, but the qualitative relationship of the efficiency to 
pitch-to-diameter or expanded area ratio in the modeling 
equations is yet to be determined. That the points of minimum 
dissipation or points of best etficiency can be determined 
from experimental data or inferred from other considerations 
is accepted, but the use of the modeling equations to estab- 
lish a ‘best’ efficiency machine given specific operating 


requirements is a topic that begs further work. 
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CHAPTER 4 PREDICTION of PROPELLER PERFORMANCE USING THE 
MODELING EQUATIONS 


4.1 Comparison with very similar propellers 


With only a few of the eighteen propellers available invest- 
igated in detail regarding the specific K(F) functions, never- 
theless it is possible to extend the results to the remainder 
of the family of propellers. As seen by a comparison of pro- 
peller numbers %Y, 15, and 18, the KC) function retains the 
same general shape over the range of P/D ratios of 1.0 to 1.4 
and EAR’s of 0.8 to 1.1. Intuitively at least, the higher EAR 
appears to smooth out the Ki and K2 functions in the dissipa- 
tion regions while slightly increasing the average magnitude 
overall. The effect of EAR is more easily seen in the raw data 
plots of Appendix B. In Figure B4¢b) for example, the higher 
EAR has its most pronounced effect at F equals +90 and -#0 
degrees. This is also shown in Figures 2.2, 2.3, 2.4, and 2.5. 
A point by point comparison will also indicate that the varia- 
tion in EAR has little effect on the thrust characteristics 
other than at those two particular points of the 4-quadrant 
curve. Seen more clearly in Figure B4¢a), it is reinforced in 
Figures 2.2 through 2.5. The expanded area ratio has as some- 
what larger effect on the Ki curve (torque), particularly in 
the normal turbine region of the 4-quadrant curve. This is not 
obvious in Figure B4(b) where the plot of the raw data would 
indicate that EAR has essentially no effect outside the prin- 


ciple dissipation regions. 
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The effect of changes in pitch-to-diameter ratio is much more 
obvious. Beginning with the raw data of Figure B3(a>, we can 
see the effect of P/D ratio throughout the analysis. An exam- 
ination of Figures 2.2 and 2.3 show this influence quite well. 
Note however, contrary to intuition, the cioser to the origin 
that a given point lies, the higher the relative magnitudes of 
the corresponding Ct and Cq coefficients. A comparison of 
Figures 2.5 and 2.6 reveal that the influence of P/D ratio is 
nearly the same for both the thrust and torque characteris- 
tics. Table 2.3 points out that the various principle frequen- 
cies of the KC PF) function change in amplitude as well as the 
average ‘dc’ term. The substantial variations in the peak 
amplitudes of the KC™ function are also evidenced in the raw 
data of Appendix B. The apparent critical points in Figure 
B3Ca> are intriguing. The increase in P/D ratio is reflected 
in the relative specific speeds of the various propellers and 
is demonstrated in the gradual counterclockwise rotation of 
the zero thrust and torque radials. Those apparent critical 
points must lie along the pivot axis of that rotation but the 
reason for the correspondence isn’t cliear. 

By examining Tables 2.2 and 2.3 and Figures 2.2 through 2.6 as 
well as the figures of Appendix B, it is possible to infer 
with very reasonable accuracy the appropriate values of the 
rating points and the fundamental values of KC® for all of 
the propellers of the data base. With that Knowledge, the 
thrust and torque plots could be generated and compared with 


Figures 2.2 and 2.3. 
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By plotting the higher harmonic amplitudes versus the zeroth 
harmonic amplitude, a comparison can be woul between three of 
the propellers examined in detail regarding the KC) function 
Figure 4.1 shows the amplitudes of the first and second 
harmonics plotted against the zeroth harmonic amplitude for 
propellers number 9, 1S, and 18. By extrapolating aiong the 
curves for a propeller whose pitch-to-diameter and expanded 
area ratios would suggest a particular zeroth harmonic 
amplitude, we can approximate the values of the Fourier series 
expansion terms for the first two harmonics. 

As an example, propeller number 1 has a pitch-to-diameter 
ratio of 0@.6 and an expanded area ratio of 0.5. Assuming the 
zeroth harmonic amplitude 1s equal to four, the first two 
harmonic amplitudes can be inferred by extending the curves to 
the appropriate value of the zeroth harmonic. Using the 
modeling equations with those amplitudes of the respective 
Fourier series expansion terms, a comparison can be made 
against the experimentally determined thrust = and torque 
curves. Figure 4.2 shows the result for the positive thrust 
curve. 

Owing to the influence of the higher harmonics near J equals 
zero, some error is expected in that region of the thrust 
curve. However, as Figure 4.2 indicates, even extrapolated 
values of the K(F) expansion terms yield quite reasonable 
results for the normal pump and turbine portions of the 
curves. The value of the zeroth harmonic (first term of the 


Fourier series expansion) has a substantial influence in’ the 
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neighborhood of J equals zero and acts with the neglected 


higher order terms to produce part of the error near the J 
equals zero values of Figure 4.2. With more study of the 
trends in the harmonic amplitudes with pitch-to-diameter ratio 
and expanded area ratio, better approximations of the thrust 
and torque curves of similar propellers will be possible. 

Another important source of error in Figure 4.2 is in the plot 
of the experimental data. The scaling which was done to 
produce the rationalized data depends on the values of n and q 
at the point of minimum dissipation. Those scale factors are 
quite sensitive to errors in the location of those points, 
sufficiently so to cause errors much like those of Figure 4.2. 
A comparison of the curves of Figure 4.2 also points out a 
weakness of the Fourier series representation as opposed to 
other possible expansion methods in forming the KCF) approxi- 
mation. The absence of the sharper peaks in the curve gener- 
ated by the modeling equations indicates that a fairly large 
number of terms of a Fourier series expansion is required to 
reproduce the experimental data in the fourth quadrant dissi- 
pation region. A piecewise continuous approximation of KCF) in 
the fourth quadrant would provide improved results for a 
simulation requiring better correlation to experimental data 


during say, transient analysis. 
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4.2 Extrapolation to other propellers 


As developed in this thesis, the modeling equations are quite 
general and although they owe their establishment to a parti- 
cular set of three bladed propellers, they are equally appli- 
cable to other types of propellers as well. Hecker [21] pre- 
sents 2@-quadrant data for a few seven bladed propellers 
including a hybrid. Plotting the experimental data in the same 
manner as with the three bladed machines, one is struck with 
the similarity of form of the curves to those presented earl- 
ier. Although not shown here, a comparison of the experimental 
data to that generated by the modeling equations using only 
approximate values for the various parameters was quite suc- 
cessful. Although the principle areas of error were in the 
dissipation regions as expected, the comparison did not indi- 
cate any serious problems with the applicabilty of the mode- 
ling equations. A more detailed comparison of the results 
obtained from the modeling equations to much different propel- 
lers is yet to be done. Considering the scarcity of complete 
4-quadrant experimental data for a large variety of propeller 
types and the limitations inherent in this form of presenta- 
tion, such a comparison wasn’t possible. However, the success 
of the equations to date would tend to enforce the argument 


for their wide ranging applicability. 
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4.3 Extrapolation to missing quadrants 


One of the advantages of the modeling equations lies in their 
ability to reconstruct essentially the entire 4-quadrant diag- 
ram using only partial information. Particularly when suffic- 
ient information is available to determine the point of mini- 
mum dissipation, a "universal* KC function can be scaled to 
fit the rationalized data and the appropriate values of Vp and 
ND at the point of minimum dissipation substituted to generate 
either the nondimensional Ct and Cq values or the dimensional 
values of thrust and torque directly. If only the first quad- 
rant point of minimum dissipation is Known, the third quadrant 
can be approximated with little loss of accuracy. This 15s 
particularly useful where first quadrant data exists but the 
model from which the data was taken is no longer available. 

Extending the concept slightly, a comparison of available data 
on the rating points of various propellers could form a_ data 
base from which the approximate rating point of a specific 
machine could be inferred. Thus, if only very limited experi- 
mental results are available, the 4-quadrant curve could be 
constructed with some confidence by comparing the modeling 


equations results with whatever information is available. 
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CHAPTER 5 SUMMARY 


5.1 An overview of the method 


The application of Eulerian turbomachine theory to the model- 
ing and simulation of pump/turbine systems has a long = and 
successful history. While the particular method used here 15 
somewhat more recent, this thesis represents an initial at- 
tempt to apply the methodology to high specific speed extended 
machines. Using very complete experimental data for a family 
of three bladed propellers with varying expanded area and 
Ppitch-to-diameter ratios as a basis, this procedure was able 
to reveal some of the important fundamental behavior under- 
lying the dissipation mechanisms of propellers in general. 
Perhaps more importantly, the model equations were able to 
Point out a common dissipation mechanism for thrust and tor- 
que, 1.e., the mechanisms which contribute to the overall 
energy losses of the system can be seen to ’reflect’ through 
the transformation of rotational mechanical to translational 
fluid energy. The plotting of the KC function for different 
propellers showed their strong correlation in form, a correla- 
tion which could in amore thorough examination be seen in the 
relative magnitudes of the components of the Fourier series 
expansions. Certainly, the zeroth component of the expansion 
represents an important parameter since a quite reasonable fit 
can be made over much of the four quadrants utilizing it 
solely in the determination of the model 4-quadrant diagram. 


This method also introduces to the simulation of propeller 
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characteristics a procedure which avoids most of the 
difficulty with singular values of the advance coefficient. It 
reduces the number of required curves to represent complete 
propeller characteristics to only two, both of which are 
reasonably well behaved and quite smooth. The model] equations 
are relatively simple and basically straightforward. Thus, 
they lend themselves to computer manipulation with Jittle 
difficulty. 

From the standpoint of the ship designer, the simulation of 
ship dynamics which include a realistic propeller in the 
initial stages of the design would be much easier and certain- 
ly more cost effective. A generic set of equations with the 
specific values of a given propeller rating point incorporated 
in them would be sufficient to provide reasonable results in 
the initial design stages which are oriented more toward near 
steady state ahead and astern operation than a thorough simu- 
lation over all quadrants. Once the geometry of the propeller 
is determined, the equations can be updated to provide in- 
creased accuracy. That accuracy depends on the number of terms 
of the expansion one is willing to incorporate into the simu- 
lation program. Thus, a scale model of the propeller is not 
necessary until the need arises for a self-propelled model 
test or a detailed investigation is required into the cavita- 


tion or vibration characteristics of the propeller. 
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39.2 The potential for further applications 


Since this thesis represents an initial examination of the 
applicability of Eulerian turbomachine theory to propellers, 
little has been addressed regarding a quantitative comparison 
of various propeller types. It was demonstrated that the 
results obtained for the three bladed propellers which formed 
the data base for this analysis could be applied to seven 
bladed propellers as well but only to the extent of a quick 
Qualitative look for comparison purposes. 

There is a great deal that can be done to enhance the useful- 
ness of the modeling equations. In particular, values for the 
KC FP) function could be tabulated based on the examination of 
avallable propeller data for various combinations of governing 
parameters. For example, the first few terms of the Fourier 
series expansions could be tabulated for given § pitch-to-dia- 
meter and expanded area ratio combinations and used in a table 
look-up scheme such that a ship designer wishing to simulate a 
given propeller would need only select a particular rating 
point and combination of P/D and EAR ratios. The computer 
program could then select the necessary coefficients and pro- 
ceed with the computations. Or, if a comparative study of the 
rating points of various well designed propellers used in ship 
Propulsion were to be conducted, a data base could be estab- 
lished such that a ship designer would need only select a 
generic type of propeller to begin a simulation. The computer 


would determine the appropriate rating point and corresponding 
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coefficients to initialize the program. 

The form of the modeling equations make them particularly 
useful for applications to controllable, reversible pitch 
propellers (CRPP). Not only are they convenient for conducting 
simulations, but they could form the necessary framework for a 
control system synthesis. Such a control system could be 
useful in automated propulsion systems particularly involving 
CRP propellers. 

There are other potential uses including applications 
involving high thrust, low speed screws like those used in 
small harbor craft or tugs. The applications to exotic 
propellers with very high sKew rates or perhaps hybrid 
propellers are also worthy of investigation. Yet another use 
involves an examination of supercavitating propellers and they 
modeling for simulation. With the promising results obtained 


here, follow-on efforts could be quite rewarding. 
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Appendix A Essentials of Eulerian Turbomachines 


In a “black-box” representation, a propeller 1iKe any turbo- 
machine can be modeled as in Figure Al where the various 
parameters are torque (tt), rate of shaft rotation (NO, fluid 
velocity (€V), and fluid flow rate (F). We define an IDEAL 
propeller as one that, in the process of power conversion from 
mechanical rotation to axial fluid flow, experiences no los- 
ses. Such a propeller may have little connection with reality, 
however the analysis of such a machine can provide qualitative 
insight into the mechanics of the energy transfer process. A 
more comfortable model would include the effects of various 
types of losses into a PRACTICAL propeller mode! diagrammed as 
in Figure Az. 


The losses can be calculated by equation Al, 
Pd = TEN —- PV CAL) 


where Pd represents the power dissipation. The amount of 
dissipation varies with the advance ratio such that two mini- 
mum will occur, one in the range of normal pump operation and 
the other approximately 180 degrees away in the four quadrant 
diagram. Using the positive semi-definite dissipation function 
described by Paynter (9), the equation for calculating the 
power dissipation from the experimental data Cin nondimension—- 


al form is found to be simply: 


$(P) = 20Cqvp - CtHp (A2) 


a? 
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The same results will be obtained using either a dimensional 
or a nondimensional formulation of the dissipation function. 
By plotting the square root of the igo eation versus the 
angle equal to the arctangent of J, the quadratic nature of 
the losses will result in a nearly linear graph such that’ the 
angle at which the minimum dissipation occurs is easily dis- 
cernable. 

Figure A3S shows the behavior of the dissipation as a function 
of the advance coefficient for propeller number 1. (CInterest- 
ingly enough, of the two minimum values of “del’, the smaller 
occurs in the region of normal turbine operation.) The exten- 
sion to a PERFECT Eulerian turbomachine comes from proper 
scaling to eliminate the appropriate amount of dissipation so 
that the expression for the general dissipation function goes 
to zero at the prescribed operating point (best efficiency or 
minimum dissipation). The PRACTICAL Eulerian turbomachine then 
supplements this dissipation function by an amount equal to 


the minimum loss. 


Per Unit variables as used throughout this paper are defined 
as the nondimensional ratio of a parameter value to some 
predetermined base value. For a pump, those base values are 
taken to be at the point of best efficency or point of minimum 
dissipation. When presenting such data in a per unit ‘form, 
obviously the reference operating point must be Known to 


reconstruct the original experiment. 
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Using per unit variables has significant advantages. First, 
they allow free manipulation of the applicable equations with- 
out regard to dimensionality. Secondly, he CO aiiee of the per 
unit variables can be made to range in the neighborhood of 
easily handled numbers such as +1 and -1 by proper selection 
of the base values. Thus the per unit values give a qualita- 
tive as well as quantitative feel for their behavior. 

Since the propeller data which provided the basis for the 
derivation of the propeller equations was already in a non- 
dimensional form such that the physical dimensions of Minio- 
vich’s propellers was no longer readily apparent, the analysis 
was done in a non-dimensional form to retain its generality. 
This may prompt one to conclude that the computational results 
are somewhat arbitrary, however one need only select Vp, WN, 
and D to derive all the appropriate values. As utilized within 


the paper, the per unit values are found by: 


P=P/Po , t=T/To , n=N/No , q=F/Fo , etc. 
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Appendix B Selected Experimentally Determined Propeller 
Data Converted to “’4-Quadrant’ Form 


The data presented by Miniovich [10] and converted to _ the 
modified advance coefficient form [C11] was converted to 4- 
quadrant form by plotting Ct (thrust coefficient) and Cq 
Ctorque coefficient) in polar form utilizing the convention 
shown in Figure Bi. This 1s possible by recognizing that at a 
constant thrust or torque Cwhich the Karman-Knapp Circle dia- 
gram utilizes), the defining relations used by Baker’ and 


Patterson shown be! ow: 


Ct = T / PD2(Up2+N2D2) 


Cq = OQ / PD2(Up24+N2D2) 


we = Arctan Vp/ND 


can be rewritten in per unit form as: 


R= 1/7 SQRTC ABSCCt) > for the thrust curve 
and 

R= 1/ SQRT¢C ABS(Cq> ») for the torque curve 
and ew = aArctan q/n for both. 


AS a result, the overlap of the data in the ranges of 


—1.5<¢J¢-0.667 and 0.667€1/3JC1.5 act as a partial check by 
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plotting on top of each other. This Kind of representation of 
the data yields in two closed curves the complete characteris- 
tics of the propeller. Wylie and Streeter £22] address this 
method as applied to solution of transient problems using a 
digital computer. An additional advantage is the avoidance of 
the problems with singularities where ever propeller speed of 
advance or shaft rpm approaches zero. See also Reference [23]. 
The curves are easily recognizible as typical for high speci- 
fic speed, axial flow pumps. See Reference [24] where Knapp 
also describes the 4-quadrant representation of pumps and 
turbines and outlines a procedure by which the representation 
could be used to deduce the transient behavior. Figure B2 is a 
sketch of the flow vs speed plane indicating the various 
regions of operation. The zero head and torque radiais are the 
asymptotes to which the thrust and torque curves approach for 
a given machine. 

Two series are shown here. Figures B3¢a) and B3¢b> are plots 
of thrust and torque for propellers 1,3, and 6. Figures B4¢(a) 
and B4¢b) are thrust and torque curves for propellers 3,? and 
iS. The plots reveal the effect of variation in pitch-to- 


diameter ratio and expanded area ratio. 
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Appendix C The Bond Graph Representation 


Paynter [91] presents the canonical bond graph representation 
for an Eulerian turbomachine as: 

T Pp 

n RG q 


It is also represented as a modulated qyrator series adjoined 


to a shunt resistor: 


Rs 





— 


_ a ee _ 
n ae Se 


In the first case, the system equations are easily found to 


be: 


J 
il 


Rg ¥ q+ Rw kn (DD 


Rg * n - Rq & q (D2) 


v0 
it 


using the pump sign convention. 





In the second case, the equations are found to be 


3 
ll 


Rg ¥ q + Rs ¥ (n-q) (D3) 


p= Rqg ¥ n - Rs & (n-q (D4) 


By converting (D3) and (D4 to canonical form, it is evident 


that Rs can be expressed as K(F) * v where 


< 
il 


SQRT ¢€ nk¥n + Qkq dD (DS) 
and 
Rg = RG + Rs & Cnkn + QkqD/VvkvV (Dé) 


A great deal of simplification is afforded with little loss in 
quality by making the simple approximation of Rg = n and Rs = 
K%¥V where K absorbs the individual character of a particular 
machine. 

AS a more widely applicable model, general dissipations could 
be adjoined to the t-junctions to account for fluid and mech- 


anical losses. i1.@e. 
Rs 


| 
ox i 


2 
ay q 


Rd Rdq 


Ce 





The equations in matrix form: 


m Rs+Rdn heh |” (D7) 


p Ro-Rs Rs-Rdq ‘la (D8) 


A more complete model would leave the dissipation elements as 
functions of the advance coefficient to absorb part of the 
individual behavior of a machine and smooth out the require- 
ments in the Rs term to fit experimental data. 

The reader interested in pursuing the bond graph representa- 
tion is directed to Reference (25] for an ‘entry level’ ap- 


proach to bond graphs and to Reference [26] for amore advan- 


ced treatment. 
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